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Conservation Genetics of an Endemic Plant, Anemoclema 
glaucifolium , in the Jinsha River Valley 


GUAN Meng-Meng'*, MA Rui’, GONG Xun!” 
(1 Key Laboratory of Economic Plant and Biotechnology, Kunming Institute of Botany, Chinese Academy of Sciences , 
SY E ) 3i 
Kunming 650201, China; 2 University of Chinese of Sciences, Beijing 100049, China) 


Abstract; Anemoclema W.T. Wang, a monotypic genus of Ranunculaceae, is endemic to the Jinsha River Valley in 
southwest China. Combining field investigation with molecular marker analysis (SNPs) , we investigated the genetic 
diversity and genetic structure of A. glaucifolium. Analysis of three chloroplast DNA (cpDNA) regions (rpsl6 in- 
tron, psbA-irnH intergenic spacer and trnC-ycf6 intergenic spacer) revealed a low level of genetic diversity within 
the species, but high divergence among populations. This genetic structure is possibly caused by a long period of 
historical geographic isolation, a relatively narrow distribution range and limited gene flow due to habitat fragmenta- 
tion. Hydropower stations are scheduled to be built in the Jinsha River drainage system in areas that overlap the 
range of A. glaucifolium and as a result of their construction the habitats of A. glaucifolium will be flooded or adverse- 
ly affected in other ways, thus threatening the survival of the species. The results of our analysis of genetic diversity 
in A. glaucifolium are of value for developing an appropriate conservation strategy for this vulnerable species. 
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Jinsha River Valley locates at the basal zone of 25 global biodiversity hotspots as featuring excep- 


the Hengduan Mountains region, which is one of the tional concentrations of endemic species ( Myers et 
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al., 2000). It is a special and vulnerable ecosystem 
characterized by aridity, high temperatures, semi- 
savanna vegetation and relatively few plant species 
(Jin et al., 1994). However, there are especially 
rich for endemic genera (e. g. , Anemoclema, Noue- 
lia, Musella , Ostryopsis, Trailliaedoxa) and species 
(e. g., Munronia delavayi, Aristolochia delavayi, 
Vitex dulouxii, Cotinus nana, Mastixia microcarpa ) 
(Jin et al., 1994). Due to natural, historical and 
anthropogenic factors, this ecosystem has become 
seriously degraded, which has resulted in soil and 
water losses and significant difficulties in vegetation 
restoration. 

Jinsha River, the upper part of the Changjiang 
River, is one of the largest rivers in southwest Chi- 
na, which originates from the Qinghai-Tibet plateau. 
A total of 14 large-scale hydropower stations are 
planning to be built along Jinsha River drainage, of 
which, some have been building, such as, Xiangjia- 
ba hydropower station, Xiluodu hydropower station 
and so on. The construction of hydroelectric power 
stations along Jinsha River will cause the rise of wa- 
ter level, so that various types of ecosystem would 
suffer serious destruction , and vertical climate would 
disappear. Moreover, many species would extinct 
owing to the submerged habitat in Jinsha River Val- 
ley. Firstly, the habitat of some species would be 
submerged, so that these species, especially endemic 
species, might extinct. For example, the construc- 
tion of the Three Gorges Dam caused the native area 
of Myricaria laxiflora to be submerged. Fortunately , 
ex situ conservation of M. laxiflora in advance pre- 
served some individuals ( Wu et al., 1998; Wang et 
al., 2003 ). 


would be destroyed after dams retained water. Verti- 


Secondly, various ecological types 


cal climate would be changed, especially the tem- 
perature and humidity of the valley. This change 
would threaten the existence of species that had a- 
dapted to dry-hot or dry-warm environment. Last, 
there would be some direct influence on species in 
construction area. The original vegetation covered 


area would be occupied by construction and living 


facilities, so that the area might lose ability to regu- 
late the ecological environment. 

Tf there is no effective measures to conduct con- 
servation strategy and population reconstruction , 
large-scale construction of hydropower stations along 
Jinsha River would cause some species to extinct in 
a short time, especially those populations close to 
the river. Studies show that before implementing 
conservation measures, genetic diversity and genetic 
structure of the species must be studied and potential 
benefits and risks of conservation strategy must be 
analyzed. Results of Marshall and Brown’ s study 
show that in order to preserve species and maintain 
normal evolution, more than 95% of the genetic di- 
versity must be preserved ( Marshall and Brown, 
1975). 

Anemoclema W.T. Wang, a monotypic genus of 
Ranunculaceae, is endemic to Jinsha River Valley in 
southwest China ( Zhang and Gong, 2002). Anemo- 
clema glaucifolium (Franch. ) W.T. Wang is peren- 
nial herb, 45-80 (-150) cm tall. The blossom is in 
July to September (Liu, 1980) (Fig. 1). A. glaucifo- 
lium grows on the rock slope of Jinsha River Valley, 
where the climate is dry and hot or warm and the alti- 
tude ranges from 1 500 m to 3 000 m. As the common 
character species of this area, it is valuable in the 
study of phylogeny and adaption to specific environ- 
ment of plants (Jin, 2002). Due to high ornamental 
value, A. glaucifolium may be cultivated as orna- 
mental or act as breeding materials after being accli- 
mated. Because A. glaucifolium is a monotypic ge- 
nus species, there are rare relevant researches except 
a few studies about molecular systematics and cytolo- 
gy (Zhang and Gong, 2002; Wang et al., 2009). 

Considering from the level of the state to deve- 
lop Jinsha River’ s water resources, it is urgent to 
protect threatened species. Detailed analysis of the 
levels and spatial distribution of genetic diversity is 
important for the development of effective conserva- 
tion strategies and management practices for endan- 
gered species (Hedrick and Miller, 1992). Therefore , 


a conservation genetics research of A. glaucifolium 
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before construction of hydropower stations on Jinsha 
River is indispensable. The research explored the 
distribution of A. glaucifolium by field investigation 
and sampling, and studied the genetic diversity and 
genetic structure through SNPs (Single Nucleotide 
Polymorphisms). Then we could elaborate genetic 
background and the distribution of genetic variation 
within and among populations, so that effective con- 
servative strategy can be made. The research can not 
only protect the endemic species A. glaucifolium , 
but also provide experience and conservative mea- 


sures to protect other biological group in this area. 





Materials and methods 
Population sampling 

In our study, leaf tissue was collected from 133 
individuals, representing 7 populations of A. glauci- 
folium from its distribution range along Jinsha River 
during July in 2008 to September in 2009. All popu- 
lations were collected in Yunnan Province, except 


for population QS, which was collected in Sichuan 


province. Information about each sampling location 
Fresh and 
healthy leaves were dried with silica gel and stored 
at 4 C. 


is presented in Fig. 2 and Table 1. 





Fig. 1 Anemoclema glaucifolium in wild. (A) The leaves of A. glaucifolium; (B) The flower of A. glaucifolium 
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Fig.2 Map showing the distribution of cpDNA haplotypes in the species. The pie sizes of sampled 


populations are proportional to their sample sizes 
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Table 1 Information of sampled populations of A. glaucifolium 
Population Locality Sample Latitude(N) Longitude( E) Altitude Haplotypes l 
Code size yS 7° /m (no. of individuals ) 
1 Qs Qiansuo, Yanyuan, Sichuan 29 28.340 101.258 2 650 Hap2 (2), Hap3(1) 
2 YS Around Yongsheng, Yunnan 10 27.374 100. 726 2150 Hapl (1), Hap2(2) 
3 SNC Shuinichang, Ninglang, Yunnan 26 27.724 101.365 2 390 Hap2 (4) 
4 LGH Luguhu, Ninglang, Yunnan 12 28. 083 101.150 2 700 Hap2 (3 ) 
5 HLT Xiangshan, Lijiang, Yunnan 24 27. 399 100.517 2 530 Hap4 (3 ) 
6 DJ Daju, Yulong, Yunnan 17 27.799 100. 484 1 860 Hap2 (3 ) 
7 ZD Hutiaoxia, Xianggelila, Yunnan 24 27.699 100.217 2 300 Hap1 (2), Hap2(1) 


DNA extraction, PCR amplification and sequencing 

Genomic DNA was extracted following the 
CTAB protocol ( Doyle and Doyle, 1987). After 
preliminary screening of a range of nucleic and or- 
ganelle DNA, we chose three cpDNA fragments 
(rpsl6, psbA-trnH, trnC-ycf6) for the full survey 
because they contained the highest number of poly- 
morphic sites. 

PCR amplification was carried out in a 20 pL 
reaction volume, containing 30 ng template DNA, 
2.0 uL 10xTaq Buffer (25 mM), 1.5-1.6 uL 
MaCl, (25 mM), 1.2-1.5 uL dNTPs (2.5 Mm), 
1.0 uL dimethyl sulfoxide (10 uM) , 0.3-0.35 uL 
each primer (10 uM), 1.5 unit of Taq polymerase 
and double-distilled water. PCR was performed in a 
T1 thermocycler ( Biometra, Göttingen, German) , 
with an initial denaturation at 94 °C for 3 min, fol- 
lowed by 30-36 cycles of denaturation at 94 °C for 
45 sec, annealing at 52-55 °C for 30 sec to 70 sec, 
extension at 72 °C for 50 sec to 1.5 min, and a final 
extension cycle of 7 min at 72 °C. All PCR products 
were purified directly using a PCR product purifica- 
tion kit (Sangon, Shanghai). Purified PCR products 
were sequenced in both directions with the same 
primers for the amplification reactions, using a 
3730xl DNA analyzer by Sangon Biotech ( Shang- 
hai) Co., Ltd. 

Sequences were aligned using Clustal X version 
1.83 (Thompson et al., 1997) and manually adjus- 
ted using BioEdit 7.0.5 (Hall, 1999). Indels were 
treated as the single mutation as substitution ( Caice- 
do and Schaal, 2004). DnaSP version 4.10 ( Rozas 


et al., 2003) was used to calculate: (1) the num- 


ber of haplotypes and variable sites; (2) nucleotide 
diversity per site (m) (Nei and Li, 1979); (3) 
haplotype diversity (H,) (Nei and Tajima, 1983). 
Within-population diversity (H,;), total diversity 
(H,), and two measures of population differentia- 
tion Gg, and Ng, were calculated using the HAP- 
LONST. Program Arlequin version 3. 1 ( Excoffier et 
al., 2005) was used to conduct an analysis of mo- 
lecular variance (AMOVA) (Excoffier et al., 1992) 
and thus to estimate genetic variations within and a- 
mong populations. Genealogical haplotype networks 
(with 95% most parsimonious connection limits ) 
were constructed using TCS version 1. 21 ( Temple- 
ton and Sing, 1993; Clement et al. , 2000). To in- 
fer possible demographic expansion of A. glaucifoli- 
um, mismatch distribution analysis, based on the 
sudden population expansion model ( Rogers and 
Harpending, 1992) using the observed number of 
differences between pairs of haplotypes, was conduc- 


ted with DnaSP. 


Results 
Nucleotide and haplotype diversity 

The aligned cpDNA rps16 data matrix was 801 
bp in length and contained two polymorphic sites 
that resulted in two haplotypes. For the psbA-irnH 
region, the aligned length was 354 bp, with three 
haplotypes derived from two polymorphic sites. Re- 
gion irnC-yc/6 was 429 bp in length, producing two 
haplotypes with only one mutation. 

The total length of the combined rps16, psbA- 
irnH , trnC-ycfd was 1 584 bp long with four substi- 
tutions (84, C/A; 525, A/C; 972, T/G; 1453, 
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C/A) and one indel (1138). A total of four haplo- 
types (Hapl-Hap4) were identified when all the se- 
quences were combined (Table 2). The nucleotide 
diversity (m) at the species level was 0.00053 and 
the total haplotype diversity (H,) was 0.325. In the 
TCS network of cpDNA haplotypes, the most widely 
distributed haplotype ( Hap2) was in the interior po- 
sition. The other three haplotypes were in the tip po- 
sition, originating from Hap2 with one to three muta- 
tional steps (Fig. 3). 
Table 2 Variable sites and distribution of four 
haplotypes of A. glaucifolium 


Variable sites 











Haplotype rps16 psbA-trnH trnC-ycfó 
84 525 972 1138 1453 
Hapl c A - c 
Hap2 c A T c 
Hap3 c A T A 
Hap4 A c G T c 


—: indicates deletion 


Fig.3 Network of the four cpDNA haplotypes of A. glaucifolium. 
Hap1 to Hap4 represent haplotypes. The size of each circle is 
proportional to the haplotype frequency. Each solid line represents 
one mutational step that interconnects two haplotypes. The small 
bold circles indicate inferred intermediate haplotypes not 


detected in this investigation 


Haplotype frequencies in each population and 
the distribution of the four haplotypes are presented 
in Table 1 and Fig. 2. Hap2 was widely distributed 
in every population except population HLT and fixed 
in population SNC, LGH and DJ. Hap1 occurred in 


population YS and ZD. Unique haploypes within 


population were detected in population QS and HLT. 
Population analysis 

On the population level, total diversity ( H, ) 
was estimated to be 0.566 and within-population di- 
versity (H,) was 0.286 (Table 3). The comparison 
between N., and Ger did not show a significant 
difference ( Ng =0. 435, Gop =0. 495; P>0. 05). 
For the total dataset, 90.31% of variation existed 
among populations, only 9.69% within populations , 


and a significant genetic differentiation was detected 


(Fg =0.90) (Table 4). 


Table 3 Genetic diversity and differentiation parameters 


based on the cpDNA fragment of A. glaucifolium 





Hs Hy Gsr N ST 
0. 286 0. 566 0.495 0.435 
(0. 1347) (0. 1527) (0.2518) (0. 2866 ) 


Standard errors are shown in parentheses 


Table 4 AMOVA of the cpDNA fragment of A. glaucifolium 





Š í anes TAR Fixation 
Source oi LE ss ariance ariation Index 
variation components /% 
CFs) 
Among a 
i 8.06 0.41 90.31 0.90 
populations 
ithi 
Wihin as’ cer 0.04 9.69 
populations 
Total 21 8.73 0.45 


d. f. Degrees of freedom, SS sum of squares, ** P<0. 001 


Mismatch distribution analysis 

The mismatch distribution based on the cpDNA 
haplotype dataset for the total samples was multimo- 
dal (Fig. 4) and inconsistent with the curve expec- 
ted for an expanding population, indicating a demo- 


graphic equilibrium (Pss»<0. 01). 


Discussion 
Genetic diversity 

In this study, four haplotypes were found by 
analyzing conjointly three chloroplast non-coding ar- 
eas (rpsl6 intron, psbA-trnH intergenic spacer and 
irnC-ycf6 intergenic spacer). The nucleotide diver- 


sity (m) is 0.00053 and the haplotype diversity 
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Fig.4 Distribution of the number of pairwise nucleotide differences 
for cpDNA haplotypes in A. glaucifolium. The dashed line shows 
observed values, whereas the solid line represents expected values 


under a model of sudden (stepwise) population expansion 


(H,) is 0.325, both are lower compared with herb 
plants in the similar distribution area (Chen et al. , 
2008 ). Genetic diversity is the consequence of long- 
term evolution. Its formation is the comprehensive 
result of sorts of factors, including current range, 
breeding system, seed dispersal mechanism and 
some historical causes like Quaternary glacial fluctu- 
ation, which can lead to repeated and drastic climate 
changes. All of these changes may cause bottleneck 
effect and founder effect, eventually, affecting dis- 
tribution and genetic diversity of plants ( Hamrick et 
al., 1992). Unfortunately, information concerning 
pollination mechanisms and breeding system for this 
species is scarce. So we deduce the low genetic di- 
versity of A. glaucifolium is the comprehensive result 
of limited individuals and some historical cause. Low 
average population genetic diversity (H, = 0. 286) 
was also detected. The reason is that there is only 
one haplotype in each four population out of the total 
seven populations, which also leads to high genetic 
differentiation among populations (G.,=0. 495, Ney 
=0.435). 
Genetic structure 

Concerning the geographic distribution of haplo- 
types, the highest frequency Hap2 was distributed in 
six populations except population HLT and was fixed 
in three populations (SNC, LGH and DJ). Accord- 


ing to coalescent theory, the haplotype in the interior 


position likely represents the ancestral haplotype, 
haplotype in tip position evolve from the ancestors 
(Wakeley, 2008). Therefore, inferred from the 
TCS network of cpDNA haplotypes (Fig. 4) , we can 
conclude that the widely distributed Hap2 is ances- 
tral haplotype and the other haplotypes ( Hapl , 
Hap3 and Hap4) all evolve from Hap2. Hapl and 
Hap3 originate from Hap2 with only one mutational 
step. Moreover, Hap! and Hap2 coexist in popula- 
tion ZD and YS, while Hap2 and Hap3 coexist in 
population QS. These further illustrate the evolution 
relationship among Hap! , Hap3 and Hap2. In addi- 
tion, the most southwestern population HLT fixes 
Hap4, which may result from geographical isolation 
or population expansion. But the mismatch distribu- 
tion (Fig. 4) showed that a recent population expan- 
sion was unlikely in this case, so the unique Hap4 
fixed in population HLT may be the consequence of 
long-term geographical isolation. 

AMOVA showed that 90.31% of genetic varia- 
tion existed among population while variation within 
population only accounted for 9.69% of the whole 
variation. In addition, we detected high level of ge- 
netic differentiation ( F',,=0.90) , much higher than 
the statistics studied by Petit (1999) , which showing 
average value of genetic markers of 97 plant ( For = 
0.70). Generally, species in fragmented habitat 
tend to have lower population genetic diversity and 
high genetic differentiation ( Young et al., 1996) , 
which is in accordance with our results. Seven popu- 
lations of A. glaucifolium in this study were distribu- 
ted in separate habitat of fragmentation. The land- 
scape characteristics of distribution area in Hengdu- 
an Mountains, where high mountains and deep val- 
leys are distributed alternately, block gene exchange 
among populations (Li and Li, 1993). Besides, 
seeds of A. glaucifolium do not have special struc- 
tures that facilitate long distance distribution or at- 
tachment to the animals, so the seed dispersal ability 
is limited. According to formula Nm = ( l-G )/ 
4G.,, we calculated the Nm of A. glaucifolium (Nm 
=0.255) , which is less than 1. Wright thought that 
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if the value of gene flow among populations is less 
than 1, limited gene flow is the main reason leading 
to the genetic differentiation of the species ( Wright, 
1931). 

Implications for conservation 

A. glaucifolium is an endemic species in Jinsha 
River Valley, which are in small numbers. The ge- 
netic diversity of A. glaucifolium is relative low with 
high genetic differentiation among populations. 
Small populations of rare and endangered species 
have higher risks of extinction than larger stable 
ones, especially when gene flow between populations 
is restricted ( Frankham et al., 2002). As small 
populations are more susceptible to the loss of genet- 
ic diversity caused by genetic drift and inbreeding, 
which reduces heterozygosity and the performance of 
various fitness related traits. 

According to field investigation, in some popu- 
lations (YS, LGH) , not only the population size is 
small, but also the number of adult individuals is 
rare, showing A. glaucifolium is in dangerous situa- 
tion of extinction. Moreover, the narrow distribution 
area and fragmented habitat result in limited gene 
flow among populations, so that genetic drift of this 
species occurred easily. Human activities also cause 
significant damages to A. glaucifolium, especially 
developing hydroelectric power in Jinshan River Val- 
ley. 1) The construction of hydroelectric power sta- 
tions in Jinsha River Valley will change the local cli- 
mate with greatly increasing the water surface area. 
On one hand, humidity will increase, which is 
harmful to A. glaucifolium. Because A. glaucifolium 
has adapted to drought environment. On the other 
hand, the daily and annual temperature ranges will 
narrow, which may alter the phenophase of A. 
glaucifolium. 2) The wider river will interfere with 
the exchange of genes between populations located 
on opposite sides of the Jinsha River. 3) The hydro- 
power station construction process produces a large 
amount of waste, which will change the soil struc- 


ture, influence the local ecological environment and 


affect local plant growth (Mu et al., 2010). 4) The 


current available farmland will be flooded, so the 
wilderness, where A. glaucifolium may grow, will be 
cultivated as farmland leading to the loss of the habi- 
tat of A. glaucifolium. Therefore, it is urgent to con- 
duct effective conservation strategy. 

Modern conservation biology aims at protecting 
not only species and habitats, but also genetic diver- 
sity of extant species for its evolutionary potentiality 
(Caughley and Gunn, 1996). Protection of genetic 
diversity needs to carry out corresponding measures 
according to genetic differences among populations. 
As to the development of protection strategy, genetic 
uniqueness is the main reference to identify prior 
population to protect ( Arthur et al., 2004). In our 
study, A. glaucifolium showed relatively low genetic 
diversity and high differentiation among populations. 
Population QS, HLT and YS contribute most to the 
total genetic diversity and allelic richness, and Hap 
4 in HLT is unique. Therefore, these three popula- 
tions should be given priority of protection. 

Taking into account current distribution and 
conservation status of A. glaucifolium, for population 
QS, HLT and YS, in situ conservation strategy 
would be of practical value. Haplotypes in popula- 
tion ZD is the same as YS. But population ZD is 
close to Jinsha River and will be affected by hydro- 
power station construction. Therefore, it is necessary 
to conduct ex situ conservation in population ZD, 
and a new population could be reestablished in the 
location with similar habitat afterwards. The seeds 
from all populations can be collected and preserved 
in Germplasm Bank of Wild Species in Southwest 
China. For the rest populations, in situ conservation 


is recommended. 
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